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Understanding many-body dynamics lies at the heart of many fundamental problems of mathe- 
matical physics. Even when one is not directly concerned with time-dependent phenomena, such 
' as in the study of equilibrium and non-equilibrium stationary states, one is essentially investigating 

, properties of the dynamics. It is also possible to approach spectral questions about the Hamiltonian 

[ of a quantum system starting from an analysis of the dynamics it generates. 

■ The equations of hydrodynamics, the Boltzmann equation, and the Gross-Pitaevskii equation, 

' are well-known examples that continue to receive a lot of attention. In these examples the aim 

is to describe the dynamics of a large number, N, of identical particles. More recently, other 
types of questions have been raised in quantum information theory, where one is less interested 
I in the large limit (although N may be large) and where one usually does not want to make 

the assumption that all degrees of freedom are identical and are subject to identical interactions. 
1^ , In either case, one is confronted with the task of taming the complexity of many-body dynamics, 

I and quantum mechanics adds another layer of complexity due to the role played by entanglement. 

Fortunately, the interactions in many physical systems are either of finite range or their strength 
decays exponentially (or with a large inverse power) in the distance between particles. In the past 
five years Lieb- Robinson bounds have been shown to be a powerful tool to turn this inherent locality 
J> . of physical systems into useful mathematical estimates. Lieb-Robinson bounds provide an estimate 

for the speed of propagation of signals (disturbances) in a spatially extended system and estimate 
the magnitude of signals propagating faster than this speed (the propagation bound is not absolute 
, since we are dealing with non-relativistic systems). 

I Systems defined on a lattice or, more generally, a metric graph, are simpler than systems of 

C3 ■ particles in the continuum because the number of degrees of freedom in a finite region of space 

can be bounded. So far, Lieb-Robinson bounds have only been obtained for such systems but we 
believe that similar results are possible for continuum systems and would also be very useful in 
applications. We will only consider lattice systems here, with the term "lattice" taken in the broad 
^ I sense of a discrete set of points. For concreteness and due to space- limitations, we will restrict 

^ ' ourselves to quantum systems. Similar considerations have been applied to classical systems such 

- - -' as anharmonic lattice oscillators \23\ [9l [37] but we will not discuss this further here. Our basic 

set-up is then as follows. 

Let r be a set equipped with a metric d. Associated to each x G F, there is a Hilbert space 
Tix and a self-adjoint operator defined on a dense domain C Hx- For example, we could 
have a harmonic oscillator at each site x. In other applications, consider the isotropic nearest- 
neighbor Heisenberg model in the absence of magnetic fields, we may have = 0. The dynamics 
is then entirely due to the interactions, which we will introduce in a moment. Thus, we have 
a multi-component quantum system defined on F, and for x, y G F, d{x, y) is interpreted as the 
distance between the subsystems located at x and y. We allow for the possibility that F is infinite, 
e.g., T = 1/^ . For any finite A C F, the Hilbert space of the associated subsystem T-L^^ and the 
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corresponding algebra of observables are given by "Ha = (S^xGA^a; -^A ~ 'S>^iT~^x), where 
B{'Hx) denotes the bounded linear operators over T-Lx- For an increasing sequence of finite subsets 
An t r, the algebra of local observables is given by the inductive limit Aioc = UnAA„. This makes 
sense because Aa„ is naturally embedded in Aa^^, for all m > nhy identifying A and A(8) 1, where 
11 is the identity operator on 'HAm\An- With this identification, for all A G ^loc we can define the 
support of A, as the smallest A such that A belongs to the subalgebra ^a- The local Hamiltonians 
for such a system are defined in terms of an interaction. An interaction is a mapping <I> from the 
set of finite subsets of T into ^loc such that for each finite X CT, ^(X)* = ^{X) G Ax- Then, a 
family of local Hamiltonians H^, A C F finite, is defined by 

(1) Ht = Y.Hx+Y.HX). 

When ^> is understood, we often suppress it in our notation. Since the sum in ([T]) above is finite, 
and each ^{X) is bounded. Ha is self-adjoint on Ha, and therefore the Heisenberg dynamics 

(2) ri'iA) = e^*^A^e-*iHA ^j. A e Aa , 
is well-defined. 

The Lieb-Robinson bounds depend on a combination of properties of F and If F is infinite, it is 
necessary to impose a condition, roughly equivalent to finite-dimensionality, as follows. We assume 
that there is a non-increasing, real- valued function F : [0, oo) (0, oo), with two properties: 

i) uniform integrablility: 

(3) =supVF(d(x,2/)) < oo, 

ii) convolution property: there exists a number C > such that for any pair x, y G F, 

(4) Yl ^(^(^' ^))Hd{z, y) < CF{d{x, y)) . 

For the case of F = W , one choice of F is given by F{r) = (1 + rY~^^. Then the convolution 
property holds with C = 2^"+^ Xlxer ^(I^D- Note that one can assume C = 1 without loss of 
generality (replace F by C'^F.) An important observation is that if there exists a function F on 
F satisfying i) and ii), then for any /i > 0, the function F^ defined by setting F^{r) = e~^'''F{r) 
also satisfies i) and ii) with \\F^\\ < \\F\\ and < C. For any ^ > 0, we denote by S/i(F) the set 
of interactions ^> for which 

x,yeX 

If $ G B^{T), then a Lieb-Robinson bound of the form 

(6) \\[nHA),B]\\ < 2\\A\\\\B\\C;' (^e'CMM _ i) J^J^W^^v)) 

xeX y&Y 

holds for all A G Ax, B G Ay, X r\Y = $, and t G M. If ^ > 0, the double sum can be bounded by 
an exponentially decaying factor of the form C\\F\\e~^''^'^'^'^\ which leads to a version of the bound 
in the familiar form: 

(7) II [t^{A),B] II < 2PII ||S||Ce-^('^(^'^)-^'*l) . 

Here, v = 2/i~^C^||<I'|l^ is the Lieb-Robinson velocity and one can take C = min(|X|, [F [) 
or in the case of interactions of finite range R, C = RC~^ m.in{\dX\, \dY\), where \dZ\ denotes 
the size of the boundary of Z. A brief discussion of the literature is in order here. The original 
bound of the form ([7]), albeit in a slightly more restricted setup, appeared in 1972 in [20] (see 
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also [U). In the following three decades, apart from a few works considering classical lattice 
oscillators [23] and a calculations for specific models [36], applications and extensions of Lieb- 
Robinson bounds received little attention. This changed in 2004 with Hastings' work [13] on the 
multi-dimensional Lieb-Schultz-Mattis Theorem [21]. In his paper, Hastings used Lieb- Robinson 
bounds in combination with his "quasi-adiabatic continuation" technique to analyze properties of 
ground states of quantum lattice systems and their excitations. For a self-contained presentation of 
quasi-adiabatic continuation see [2]. This technique and Lieb- Robinson bounds were subsequently 
extended and used in new applications by Hastings and collaborators [HI [HI [5l [HI [T71 [151 O IS] 
and other authors [28l[25l|29l[8l[32l[l0l[IIl[30llMl[l2l[Ml[3a[^ 

The new applications motivated extensions of the Lieb-Robinson bounds to more general systems 
as well as improvements on the basic estimates. We will mention some of these newer results below 
but space limitations will not allow us to explicitly mention all the results that have appeared in 
the past few years. First, however, we clarify the relation between bounds on commutators and the 
support of observables. 

Let Til and 7^2 be Hilbert spaces and suppose there is an A G B{'Hi "^2) such that for all 
B G B{n2) we have ||[^, l(g)B]\\< €\\B\\. Then, there exists A, G BCHi) such that \\A - A^\\ < e. 
In fact, if 7^2 is finite-dimensional can be obtained using the partial trace of A over 7^2) but 
the result also holds for bounded observables on infinite-dimensional Hilbert spaces [2]. For fixed 
A G Ax, one can apply this result with t{^{A) taking the role of A. By letting Tii = Tixr, where 
Xr = {y ^ A, d{X, y) < v\t\ + r} and H2 = 'H\\XrJ see that ([7]) gives a bound for the error one 
makes if one replaces t^^{A) by its best approximation with support contained in X^. Therefore, 
it is not surprising that Lieb-Robinson bounds can be used to prove the existence of the dynamics 
in the thermodynamic limit [3]. While Lieb-Robinson bounds do not always give the best possible 
estimates, in some cases they do, and moreover, they have been used to prove the only known 
results for anharmonic lattices [271 [1] . 

Another fundamental result that directly relies on the locality properties expressed by Lieb- 
Robinson bounds is the Exponential Clustering Theorem. It states that Hamiltonians with a 
non-vanishing spectral gap have ground states with exponentially decaying spatial correlations. 
Although the corresponding result in relativistic quantum field theory has been known for a long 
time, the non-relativistic version, with the Lieb-Robinson velocity playing the role of the speed of 
light, was only proved in 2006 [281 116] . The correlation length ^ satisfies ^ < 2v/'y + l//i, where v 
is the Lieb-Robinson velocity, 7 is the spectral gap above the ground state, and fi is the parameter 
measuring the exponential rate of decay of the interaction. When one uses a more recent version 
of the Lieb-Robinson bounds [30], the exponential clustering for gapped ground states uj of one- 
dimensional systems with short-range interactions, can be shown to hold in the following strong 
form: 

\uj{AB) -uj{A)uj{B)\ < c\\A\\ \\B\\ex.p {-d{X,Y)/C), 

where X and Y are the supports of A and B, respectively. Note that only the distance between 
the supports and not their size appears in the estimate, a feature that was exploited by Matsui in 
his investigation of the split property for quantum spin chains [24] . 

The Area Law, the conjecture that the entropy of the restrictions of gapped ground states of 
quantum lattice models to a finite volume A grows no faster than a quantity proportional to the 
surface area of A, has been proved by Hastings for one-dimensional systems in [13]. In his paper 
Lieb-Robinson bounds are used to derive an approximate factorization property of the density 
matrices of gapped ground states. Such a result can be generalized to higher dimensions |12j . but 
it does not, by itself, suffice to prove the Area Law in this context. This issue remains a topic of 
active investigation (see, e.g., [3]). 

The estimate of the correlation length in terms of the gap plays an important role in several 
other applications of Lieb-Robinson bounds, including the multi-dimensional Lieb-Schultz-Mattis 
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(LSM) theorem. A precise statement of the general multi-dimensional LSM theorem would be 
too long to fit in the space alotted, but since it was one of the first non-trivial applications of 
Lieb-Robinson bounds, it deserves to be discussed here. For concreteness consider the spin-5 
nearest-neighbor isotropic quantum Heisenberg antiferromagnet with S = 1/2, 3/2, 5/2, . . ., defined 
on a finite subset of Z'^ of the form [1, 2L] x [2L + 1]*^"^, and with periodic boundary conditions in 
the first coordinate. It is known that this model has a unique ground state [22]. The LSM theorem 
[13} [29] then provides a bound for the energy of the first excited state: 

El — Eo < C ^ , 

where C is a constant of order 1, only depending on the dimension and the coupling constant. The 
proof is by showing that a gap larger than the bound claimed above would allow one to construct 
variational states of lower energy, which is a contradiction. The properties of the variational states, 
the estimate of their energy and the proof of orthogonality to the ground state, rely on Lieb- 
Robinson bounds in an essential way. See |30[ Section 5] for a more detailed outline of the complete 
proof. Here we just mention that Lieb-Robinson bounds allow one to show that, as long as the 
spectral gap is not too small, local perturbations added to a Hamiltonian modify the ground state 
only in a neighborhood of the perturbation. See [2] for a general proof of this property. 

In an impressive application of the adiabatic continuation technique and Lieb-Robinson bounds 
Hastings and Michalakis proved the quantization of the Hall conductance for a general class of 
models of interacting fermions on a lattice in |17] . An extension to the fractional quantum Hall 
effect is also discussed in this work. 

Bravyi, Hastings, and Michalakis [3 [6] have shown that topological order in the ground states 
of a class of Hamiltonians that are the sum of commuting short-range terms, such as Kitaev's toric 
code model [19], is stable under arbitrary sufficiently small short-range perturbations. This can be 
regarded as another instance where Lieb-Robinson bounds are used to show that local perturbations 
have only local effects, and therefore cannot destroy a global property such as topological order. 

In this short review we have only considered Hamiltonian quantum dynamics. We just note that 
Lieb-Robinson bounds have also been derived for irreversible dynamics described by semigroups 
with a generator of Lindblad form |33j . 

The complexity of quantum dynamics and its tendency to create entangled states are a bar- 
rier to our intuitive understanding of many of the most interesting physical phenomena. Good 
mathematical results that elucidate the structure behind this complexity are essential to aid our 
understanding. We hope to have convinced the reader that Lieb-Robinson bounds are a good 
example of this kind of mathematical result. 
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